Electronic transport in graphene: A semi-classical approach including midgap states 
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Using the semi-classical Boltzmann theory, we calculate the conductivity as function of the carrier 
density. As usually, we include the scattering from charged impurities, but conclude that the 
estimated impurity density is too low in order to explain the experimentally observed mobilities. We 
thus propose an additional scattering mechanism involving midgap states which leads to a similar 
k-dependence of the relaxation time as charged impurities. The new scattering mechanism can 
account for the experimental findings such as the sublinear behavior of the conductivity versus gate 
voltage and the increase of the minimal conductivity for clean samples. We also discuss temperature 
dependent scattering due to acoustic phonons. 

PACS numbers: 72.10.-d, 72.15.Jf, 72.15.Lh, 65.40.-b 



I. INTRODUCTION 



in an insulator such as SiC>2. 16 
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Electronic transport through two-dimensional 
graphene sheets has attracted incessant attention 
ever since the first experimental measurements on a 
Hall geometry were performed by Geim and coworkers 
about three years ago 1 -. This is due to some spectacular 
findings like the universal minimal conductivity at the 
Dirac point and the high mobility of the samples which 
is basically independent of doping and temperature^^ 
Different experimental setups, i.e., electrical field 1 versus 
chemical doping^, also give rise to different theoretical 
models. For recent qualitative reviews on both the 
experimental and theoretical status of the field, sec Ref. 

Mm. 

So far, the scattering mechanism which determines 
the transport properties has not unambiguously been 
identified, but there is strong evidence that long-range 
Coulomb scatterers can account for many of the exper- 
imental findings*^ ' 10 ! 11 It was shown within the Boltz- 
mann formalism that the conductivity scales linear with 
the carrier density if one assumes charged impurities in 
the SiOvsubstrate close to the graphene-sheet^ This 
semi-classical approach was also applied to systems close 
to the Dirac point and in the presence of adsorbed 
molecules 

Recently, it was claimed that a Boltzmann theory with 
long-range Coulomb scatterers can account for all exper- 
imental findings if one renormalizes the carrier density 
close to the Dirac point due to potential fluctuations . 13 i 14 
The theory predicts a non-universal behavior of the min- 
imal conductivity at the Dirac point which nevertheless 
coincides with the experimentally observed value of 4e 2 //i 
for "dirty" samples. This is in contrast to numerical 
studies based on the Kubo-formalism by Nomura and 
MacDonalcU^ who show that the conductivity is a func- 
tion of n/rii with the impurity density, thus finding 
universal behavior when the carrier density n goes to 
zero. The main criticism of Refj 1 ^ is the high density of 
charged impurities 7ii>10 12 cm -2 needed to match the ex- 
perimentally observed mobilities, not likely to be present 



In this article, we propose a new scattering mechanism 
originating from midgap states which may be formed due 
to vacancies, cracks, boundaries or impurities in the sub- 
strate with a high potential difference with respect to 
the graphene sheetj 17 ' 18 They also occur in corrugated 
graphene* 1 ^ The phase shift resulting from these types 
of disorder must approach zero for wave vectors close to 
the Dirac point. In contrast to the phase shift due to 
a short-range contact potential, this behavior is not lin- 
ear but logarithmic ! 20 ' 21 The resulting scattering time is, 
therefore, proportional to k up to logarithmic corrections. 
It is interesting to note that this behavior is also found 
for a two-dimensional ( "non-relativistic" ) electron gas 2 ^ 
and in corrugated graphene 22 where the focus was laid 
on the resulting random gauge field. 

Within the Boltzmann approach, the new mechanism 
can account for: a) quasi-universal minimal conductivity 
for dirty samples; b) a higher minimal conductivity for 
cleaner samples; c) sub-linear behavior of the conductiv- 
ity as function of the gate voltage. We further obtain 
realistic values for the mobility assuming an equal con- 
centration for Coulomb scatterers and vacancies of order 
m ~ 10 10 - 10 n cm- 2 . 

The paper is organized as follows. In Sec. UH we 
will first introduce the Boltzmann approach and com- 
ment on its applicability to graphene, i.e., to chiral Dirac 
Fermions. In Sec. IIII1 we discuss the density of states 
in the presence of midgap states, needed to estimate the 
transport properties close to the Dirac point. In section 
IIV1 we calculate the relaxation time and electrical con- 
ductivity for the various scattering mechanism includ- 
ing acoustical phonons. In section [V] we discuss the 
ac-conductivity, the thermal conductivity and the ther- 
mopower for the new scattering mechanism including 
midgap states and present numerical results in section 
IVI1 We close with conclusions and remarks. 
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II. BOLTZMANN EQUATION 
A. Collision- free Boltzmann equation 

We start by showing that the Boltzmann equation de- 
scription leads to the same plasmon spectrum as the 
more used many-body methods. 23 ' 24 This shows that a 
semi-classical approach for the transport properties of 
graphene is accurate. 

The Boltzmann equation is described in terms of the 
electronic distribution function f k . Within this semi- 
classical approach, f k depends on space r and time t, 
i.e, 



fk = fk(t)(r,t) 



(1) 



Looking at time scales shorter than the lifetime of the 
quasi-particles, the number of quasi-particles in the state 
k is conserved. Via the continuity equation f k + V r • 
jk = with j k = v k f k denoting the particle current, one 
arrives at the collision-free Boltzmann equation. With 
k = eV r c/? where (p is the scalar potential of the internal 
electrical field, this reads in Fourier space as 26 

Q 1 

(-iu + iq ■ v k )f k (q,u)) = ieq ■ v k (--—-)ip(q,ui) . (2) 

oe k 

To investigate screening properties, an external poten- 
tial ip ext (q,oj) is assumed. To linear order in the total 
potential tp, the induced density is then given by 



p md (q,Lu) 



E 



q ■ v k 



A ijj - q ■ v k 



(3) 



where A is the area of the graphene sheet and spin and 
valley degeneracies have been included. 

The induced potential is obtained by (— e)tp md {q, u) — 
V q p md (q,oj) with V q = the two-dimensional Fourier 
transform of the Coulomb potential. For the dielectric 
function e(q,u>) = ip ext (q,u>)/<p(q,uj), one then obtains 
in the long-wavelength limit vpq <C to the following ex- 
pression: 



e(q,u) s» 1 



"« 4 B*--) 2 (- a/2 



2 A 



de k 



= 1 - 



V q q 2 e F 



(4) 



Plasmon excitations are given by e(q, u) = which leads 
to the plasmon dispersion 



27re 



eFq 



(5) 



This relation including the prefactor is also obtained from 
a standard tight-binding model of graphene where the 
dielectric function is calculated within the random-phase 



approximation 



23,24 



Our subsequent results should thus 



be valid even close to the neutrality point as long as 
kp£ ^> 1 (£ the mean free path), i.e., the chirality of 
the Dirac fermions only enters in the expression for the 
transition rate (see Eq. (|2~lj) ). For a quantitative analysis 
starting from a two-band model, see Refi2£. 



B. Collision term 

We now include the possibility of changing the quan- 
tum state k by introducing a collision term which is usu- 
ally facilitated by the relaxation-time approximation: 27 



dfk 

at 



<jk 



(6) 



where fh~ f k = 9k- 

Applying an electric field E to the sample, the solution 
of the linearized Boltzmann equation then reads 



9k 



df(e k ) 
de k 



er k v k ■ E . 



and the electric current reads 
. 4 



-7 E eV k9k ■ 
k 



(7) 



(8) 



Since at low temperature the following relation 
—f°(tk)/de k — > 5(vph(k — kp)) holds, we obtain for the 
conductivity with the Fermi velocity vf the well-known 
formula^ 



-p(E F )T kf 



(9) 



In the following, we will give expressions for the density of 
states p{E) and the relaxation time T k - We then discuss 
the electrical conductivity in the low- and high-density 
limit. 



III. DENSITY OF STATES 

The density of states per unit area of clean graphene 
is given by 



Po(E) 



2\E\ 



(10) 



where spin and valley degeneracies have been included. 

Due to potential disorder this linear behavior becomes 
sub-linear—, though the density of states at the Dirac- 
point is still zero. More important are local defects in 
form of vacancies which were first discussed in Ref. l30l . 
Within a self-consistent Born approximation (CPA), it 
was shown that the relaxation time depends linearly on 
the mean free path, i.e., r ~ £/vp- For the mean free 
path we have £ ~ 1 /■y/rii, where rij stands for the impurity 
density due to vacancies, cracks, etc. . In order to obtain 
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an universal minimal conductivity, the density of states 
close to the Dirac point must be given by p <~ \J(lhvp). 

This behavior is also obtained from a phenomenolog- 
ical approximation. For this, it is important to note 
that vacancies gives rise to localized states which decay 
algebraically^ These states hybridize due to the over- 
lap with localized states of other vacancies. The energy 
scale is given by the mean distance between vacancies 
and approximated by the gain of energy due to the new 
boundary conditions. 

This energy scale is approximated as follows. The lin- 
earized tight-binding Hamiltonian for a graphene sheet 
with circular symmetry is given by 



H, = 







a is<f> 



s*(_ isdr _l d(j> ) 



(~isd r + ±84) 




(11) 



where s = ± denotes the two valley. At non-zero energy, 
the general solution is given by the Bessel and Hankel 
functions. Considering only one valley 8=1, the general 
wave-function in graphene at low energies is thus given 
by 



Jo(kR) 
-U l {kR)e i < l > 



B 



Y (kR) 
-iY 1 (kR)e l ' t > 



(12) 



A simple model for vacancies of radius i?o and localized 
on the radius Ri now assumes that the first component 
becomes zero at the inner boundary Rq and the second 
component at the outer boundary Ri , thus assuming zig- 
zag edges on different sublattices. This leads to the fol- 
lowing quantization condition for k: 

Jo(kRo)Yi(kRi) - Yo(kRo)Ji{kRt) = (13) 

For kRi -C 1, the lowest momentum is then given by 

11 

(14) 



Ri y/\HkRo)\ 

which defines the width of the localized band 



E 



loc 



fliVIKRo/fli)l 



(15) 



The density of states at zero energy is thus approximated 
by 



p(0) = 



V\HR /R 1 

HvpRi 



(16) 



Since i?i is related to the average distance between the 
vacancies, we have p(0) ~ (ru\ In nil) 1 / 2 . Notice that 
the CPA calculation of Ref. [30 does not capture the 
logarithmic correction. We will therefore approximate 
the density of states as 



(17) 



p (E)=a^g c (E loc -\E\) + ^p^ 
nvp Tr(HvF) 



where a « 1/2 is a dimensionless constant and g c (E) 
is a cutoff function, e.g. g c (E) — 0(E). We note again 
that we assume coherent impurity scattering at work, i.e., 
p(E — > 0) ~ y/nl- For clean samples, we would expect 
the standard scaling behavior, i.e., p(E — > 0) ~ rij. 

The density of states of Eq. ((TTJ) characterizes two 
regimes. For high carrier density Ep > E\ oc , the conduc- 
tivity of Eq. © reads 



2e 5 



VFk F Tk f 



(18) 



Close to the Dirac point Ef < E\ oc , we obtain the mini- 
mal conductivity 



27re 2 

Cmin = —r—VF{a^Jni)T kf 



(19) 



Notice that we obtain the same formula as for high elec- 
tronic densities by introducing a minimal Fermi wave- 
n7. The minimal Fermi wave- vector with 



vector fc m i n 

kpl ^ 1 can be related to self-doping effects induced 
by the very same mechanism which is invoking midgap 
states^ and we will show in the following that a crossover 
from a linear to constant behavior of the conductivity 
versus gate voltage takes place. 

If one believes that close to the neutrality point the 
system also behaves in a diffusive way, i.e., that the ex- 
perimentally observed negatively and positively charged 
puddles form a macroscopic network^ then even for low 
impurity densities rii ~ 10 10 cm -2 with kpi ~ 1 one can 
use our estimates (up to a constant of order 1). Still, 
we cannot rule out the existence of another regime where 
the Boltzmann approach is invalid and e.g. percolation 
models are at work^ 2 - 

Having established the typical behavior of the density 
of states at the Dirac point due to midgap states, we will 
discuss several scattering mechanisms. 



IV. RELAXATION TIME AND DC 
CONDUCTIVITY 

The collision rate I/77- due to impurity scattering is 
usually given by 2 ^ 



1 

Tfe 



= JVi £ r(fc, *')(!- cos fc ,tf 



(20) 



where N% is the number of impurities and the transition 
rate from the quantum state k to k' is approximated by 
Fermi's Golden rule 



2tt 



T(k,k') = — \{k\V scat tW)rS(E k -E k ,) . (21) 

It is only in the scattering matrix r(fc, k'), where the chi- 
rality of the Dirac fermions enters within the Boltzmann 
formalism. If the scattering potential does not break the 
sub-lattice symmetry, this will only lead to a numerical 
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factor. With the Fourier transform of the scattering po- 
tential V sca tt(q), the collision rate can then be written 
as 



p(E F ) / d9\V sc Mq)\ 2 (l - cos 2 9) , (22) 



h n 

T kF 

where ri| catt is the impurities density of the scattering 
potential and q = 2kp sin(0/2). Notice that the argu- 
ment of the integral vanishes for both 9 — and 9 = n, 
a situation that does not occur in normal metals. 

The effect of vacancies or local impurities with a high 
potential difference with respect to the graphene layer 
cannot be treated with the above formulas since they do 
not capture the change in phase space due to midgap 
states. We thus determine the relaxation time via the 
phase shift induced by the scattering center. Assuming 
elastic scattering and only considering s-wave scattering, 
the transition-rate is then expressed as2£ 



h 

Tk 



8m 



np(E k 



sin 2 (4) 



(23) 



where 8k is the phase-shift of the s-wave channel. 

In the following, we will consider various scatter- 
ing mechanisms, i.e., we will discuss the effect on the 
electronic conductivity due to a) local substitutions 
(short-range "contact" potential); b) charged impurities 
in the SiC-2 substrate (long-range (screened) Coulomb- 
potential); c) acoustic phonons, where Eqs. (f2T)|) and 
([2"Tj) have to slightly be modified. In Subsection IIV Dl we 
introduce the new scattering mechanism due to midgap 
states. 

Due to the unusually high energies of optical phonons 
of the order of 0.1-0.2eV in graphene-related materials, 
optical phonons cannot be treated within the Boltzmann- 
formalism since they induce interband-transitions for 



usual densities n<5 x 10 12 cm 



For a discussion on 



transversal optical phonons in graphene sheets within the 
Kubo-formalism, see Refi^2. 



A. Contact potential 

We will first discuss the scattering behavior from 
V^catt(f) = v o<5(l r D- This yields a relaxation time 



l 



Tfc 



4h 2 v F 1 



-.Contact 



nv^ p(Ek) n 



Co7itact„,2 



(24) 



where 



-.Contact 



is the impurity concentration and the 



right hand side resembles the high carrier density limit. 

Eq. (f2"4")) can also be obtained from calculating the 
phase shifts. From Refj2£ we obtain Sk = VQk/(Ahvp) in 
the limit of small k and by expanding Eq. (|23[) up to 
linear order in Sk, we obtain the above result. 

The conductivity does not depend on doping^ and we 
obtain 



Contact 



„ Contact 



{hv F f 



^Contact „,2 



(25) 



Eq. (j25|) does not depend on an energy scale nor does it 
lead to a significant contribution for the total conductiv- 
ity. This is more generally known as the Klein-paradox. 21 



B. Long-range Coulomb potential 

Let us now discuss the influence of the long-range 
Coulomb potential. Charged impurities reside in the iso- 
lating Si02-layer and are screened by the conduction elec- 
trons of the graphene sheet. This yields for the potential 



m momentum spao 



,35 



2e e q 2e e q 



(26) 



where p md (q) is the induced charge density, e the per- 
meability of the substrate, and z c denotes the shortest 
distance of the external charged impurity to the two- 
dimensional graphene sheet. 

Since we are employing a semi-classical approach, it 
is consistent to approximate the induced charge density 
within the Thomas-Fermi (TF) approach: 

P md {r) = -e- - Mr)) - /(<*)) 



eMrAj2(-7^) = - e2( P(r)p(eF) (27) 



A 



The last equality follows by assuming a Fermi liquid char- 
acterized by a sharp Fermi "surface" . For a discussion 
on non-linear screening, see Ref. l36l . 

The TF approach thus gives the following form for the 
screened potential inside the graphene sheet: 



v{q) 



Ze e- q \ z "\ 
2e e q + j 



(28) 



with 7 = p(Ep)e 2 /2eoe and the density of states given 
by Eq. (TJU). 

At large doping we have 7 = ^kp and from Eq. (|22"|) 
with V S catt(q) = e<p(q) and z c ~ we obtain 



H 2 VFkF 



, with uq 



4e e(l + 7) 



(29) 



where nf is the density of charged impurities in the sam- 
ple. This leads to the following conductivity: 



coulomb _ 2e2 {hv F k F ) 



(30) 



Eqs. and fl3DJl are slightly modified for z c > 03 
At zero doping, we obtain the minimal conductivity 



Coulomb 



4e 2 rii 2 

-r — c 2") 

n ny 



(31) 
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We thus find universal behavior at low doping if rii « nf . 
For a — 1/2, we obtain the the experimentally observed 
value of a = 4e 2 jh. 

Now we want to determine the numerical values of the 
relaxation times due to charged impurities and later com- 
pare them to the ones due to acoustic phonon scattering. 
Let us assume that the electronic density in the graphene 
sample, induced by the gate voltage, has the typical value 
(gate voltage of 100 V)i 

n = 7.2 x 10 12 cm" 2 . (32) 

The Fermi momentum is given by n = k 2 F /ir, where con- 
tributions from both spins and both Dirac points were 
included; this leads to a value for k F given by 

k F = 4.8 x 10 8 m" 1 . (33) 

With e 2 /47re = 14.4eVA and hv F = 3at/2 = 5.75eVA 
where a — 1.42A and t = 2.7eV, we find for large densi- 
ties 7 = 10k F /e. For Z — 1 and z c ss we set e = 2.4 
which is the average of the dielectric constant of SiC>2 
and vacuum, since the graphene-layer is sandwiched by 
these two layers. For large carrier densities, we thus have 
7 = 7&p with 7 « 4.2. 

Using Eq. (|29|) and the above value for k F we obtain 

r kF ~ 6 x HT 17 (nf y 1 s, (34) 

where nf is the concentration of impurities per unit cell. 
Since it has been experimentally determined that the 
mean free path £ of electrons in graphene can be as large 
as 0.3 fim. the experimental relaxation time is seen to be 
of the order of 



— ~ 3 x 10 

v F 



-13 



(35) 



This last result implies that the concentration of impuri- 
ties per unit cell has a value of the order of 



2 x 10" 



(36) 



if one only assumes the Coulomb scattering mechanism. 
This corresponded to a density of nf = 4 x 10 11 cm -2 
and a mobility of 11100 cm 2 /Vs. But note that graphene 
normally exhibits lower mobilities for which a higher im- 
purity density is necessary^ Also finite z c and a higher 
dielectric constant, e.g. e = 3.9 which holds for SiC>2, 
lead to a larger relaxation time and thus imply higher im- 
purity concentrations to match the experimental results. 
In summary, we strongly believe that the Coulomb scat- 
tering mechanism is only capable to explain a marginal 
fraction of the observed data. 



where the transition rate T(k, k') is given by 

r(fe, k') = ^\H k , k \ 2 S(v F hk' - v F hk - Two) , (38) 
n 

where v F hk is the dispersion of Dirac fermions in 
graphene, hut the phonon energy and H^ k is defined 
as 



Hi 



k'.k 



dr^l,(r)Us{r)ipk{r) , 



(39) 



with Us(t) the scattering potential and ipk(r) is the elec- 
tronic wave function of a clean graphene sheet. 

If the potential is due to the propagation of phonons 
it has the forrn^ 



U 3 = K„A q e l{qr - ut) 



when 



,38 



\K„ 



n 



L <3 



2pAu> q 



N(uj q 



1 



k B T 



e hu q /(k B T) „ l fjjj 



(40) 

(41) 
(42) 

(43) 



where p is the density of graphene and Da is the electron 
acoustic deformation potential, estimated to be of the 
order of 3t, where t is the first neighbor hopping matrix 
in graphene of the order of 3 eVi^ A similar estimate 
for the deformation potential is obtained by relating the 
bond-length with the hopping amplitude^. 
The matrix element Hk'.k is easily obtained as 



H k ,, k = cos(6/2)K q A q 6 

k-\-q,k' & 



(44) 



Using Eq. JH]) in Eq. l(38|l. the transition rate T(k,k') 
reads 



r(fe,fe') = -|*g 2 |A,| 2 (i + cos0)<5 fe+< ^ 

x S(v F hk' — v F hk — huj q ) . 



(45) 



The form of T(fc, k 1 ) represents the absorption of a 
phonon of momentum q and energy hw q . Since we want 
to describe the absorption of acoustic phonons, we write 
the dispersion uj q as 



u q = v s q , 



(46) 



where vs is the sound velocity. The conservation of mo- 
mentum, k + q = k' , leads to 



C. Phonons 



The relaxation time Tk for phonon scattering is defined 



— = Vr(fc,fc')(l -cos0), 



(37) 



Vfc 7 



2fc'fccos6». 



(47) 



which allows us to write the product of the Kronecker 
symbol and the Dirac delta function as a single delta 
function reading 



S(v F hk' - v F hk - TivsVk' 2 + k 2 - 2k'kcos6). (48) 
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Since vp ^> v s the argument of the delta function 1(48]) 
can be approximated by vpTik! — vphk, i.e., the absorp- 
tion of acoustic phonons can be seen as a quasi-elastic 
scattering process. 

The final result for the scattering time is 



8h 2 pVgVF 1 

D\k B T k 



(49) 



which is formally similar to the relaxation time produced 
by a contact potential, Eq. (|2~4"j). except for the temper- 
ature dependence which is absent in the latter case. 

Let us now concentrate on the numerical values of the 
relaxation time due to phonons. The phonon spectrum 
of graphene has two acoustic branches named LA, with 
a velocity of 7.33 x 10 3 m/s, and TA, with a velocity of 
2.82x 10 3 m/s4^ At the temperature of 1 K the relaxation 
times of these two modes are 



TLA 



2.7 x 1(T 10 s, 



T TA ^ 4.0 X ict 11 S. 



(50) 
(51) 



which are clearly much larger than the estimated value 
of 3 x 10~ 13 s. Therefore at this low temperatures the 
scattering is mainly dominated by impurities. At tem- 
peratures around 100 K the scattering life time dimin- 
ishes by a factor of 100 leading to values comparable to 
that obtained from the scattering from charged impu- 
rities. Considering only charged impurities as scatter- 
ing mechanism, the effect of phonons on the transport 
properties of graphene must be taken into account if the 
calculations of the transport coefficients are extended to 
temperatures of the order or above 100 K. 

Let us finally discuss the effect of temperature on the 
need of taking into account thermal excitations of the 
valence band into the conduction band. For the value of 
the Fermi momentum given above the Fermi energy has a 
value of hvpkp = 0.3 eV. This energy value corresponds 
to a temperature of the order of 3600 K. Therefore as 
long as the temperature is much smaller than this value 
the valence band can be considered as inert and therefore 
we can perform the calculations by neglecting the effect 
of the valence band altogether. 



In order to incorporate the effect of midgap states in 
the calculation of the relaxation time we depart from Eq. 
(f2"3"| . Scattering from vacancies leads to the following 
phase-shift)^ 

1 



C 

k ~ ~2\n(kR ) 
This means that for fci?o <C 1 

h P {E k ) 



n. 



-(InkRo) 



(52) 



(53) 



2irn,i 

For large carrier densities, p(Ek) ~ fc, and apart from 
the logarithmic correction, this is the scattering behav- 
ior coming from long-range Coulomb potentialsi&^i 3 ^ Ex- 
plicitly, one obtains 

r k = ^^(\nkR ) 2 . (54) 

7T VpTli 

The logarithmic correction leads to a sub-linear density 
dependence of the conductivity 



Vacancies 



^-k 2 F (\nk F R ) 2 

h 7T 



(55) 



We note that the same behavior is obtained if one in- 
cludes wiggles and thus a random magnetic field in the 
graphene sheet. = The relation between midgap states 
and corrugated graphene was investigated in Ref.— . 

1/2 

For zero carrier density, we obtain with fc m j n ~ n i the 
following result for the minimal conductivity (n,; = fi{j A c 



and i?o ~ A c ): 



Vacancies 



— (a | In ni|) 2 



(56) 



Notice that there is no "linear" dependence of the im- 
purity density n,j. Having a typical impurity density per 
unit cell of fij = 10 which matches well the experimen- 
tally observed mobility, the logarithmic correction can 
be approximated as |lnni| ~ 8. For a — 1/4, we ob- 
tain the experimentally observed minimal conductivity 

CTmin = 4e 2 //l. 

For cleaner samples the logarithmic correction has to 
be taken into account. This is in accordance to ex- 
perimental findings which show higher conductivity for 
cleaner samplesj 14 i 42 



D. Vacancies 

Vacancies, cracks or boundaries in the graphene sheet 
give rise to bound states at the Dirac point, so-called 
midgap states. This is also true for corrugated graphene. 
There is thus inherent disorder in the system that has to 
be treated adequately. The influence of midgap states is 
not captured in Eqs. (|2H)) and (f2"Tj) , where the reference 
point is given by the unperturbed system described by 
plane (spinor) waves, i.e, ipk{r) = (x\k) is the electronic 
wave function of a clean graphene sheet. This is also the 
reason why the Klein paradox is not at work which would 
render scattering from local impurities (i.e., vacancies) 
irrelevant. 



V. OTHER TRANSPORT QUANTITIES 

In Refill, predictions were made on how the ac con- 
ductivity, the thermal conductivity and the thermopower 
depend on the carrier density if one assumes the scatter- 
ing behavior ~ k. The scattering mechanism from 
vacancies, cracks or boundaries yields a different scatter- 
ing behavior, i.e., ~ fc(lnfci?o) 2 . In the following we 
give the density dependence in terms of the Fermi en- 
ergy Ep — hvpkp for the above quantities assuming this 
new scattering mechanism at work (n = E F /ir) . Measur- 
ing these quantities may then disclose which scattering 
mechanism dominates. 
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A. The optical conductivity 

Here we want to obtain the electronic density depen- 
dence of the optical conductivity of a doped graphene 
plane. Since the Boltzmann approach does not include 
inter-band transitions, the expressions obtained below 
are only valid as long as fkj < Ep with Ep — hvpkp 
the Fermi energy, where the above mentioned transitions 
are blocked by the Pauli principle. 

Our aim is to obtain the response of the electronic 
system to an external electric field of the form 

E = E e l(q - r -^ . (57) 

The Boltzmann equation has, for this problem, the form 

df°(e k ) g k dg k , 

s ev k -E = h v k ■ V r g k + -— . (58) 

oe k Tfc at 

The solution of the linearized Boltzmann equation (|58|) 
is well known 2 ?, reading 



oe k 



with 



er k v k ■ E Q 



(59) 



(60) 



The Fourier component J(uj, q) of the current is given by 

df(e k ) 



J{uj,q) = J d 2 kev k <S> q (uj,k) (-- 



de k 



, (61) 



leading in the long-wavelength limit to an optical con- 
ductivity of the form 



h u, 



~ n2 



1 



1 + i^QnEp/vo) 

^0 

+ (^{\nE F /v Q y 



(62) 

In the above equation, we defined the two energy scales 
Uq = TT 2 niH 2 vp and vo = Hvp/Rq. What should be 



E 2 F /n. 



stressed about Eq. (|62|) is its density dependence n 

4/ 



B. Thermal conductivity and thermopower 

In the presence of a temperature gradient in the sam- 
ple, the linearized Boltzmann equation has the form 



df°(e k ) 
de k 



-v k 



£fc — Ep 
T 



V r T + eE 



bsv. 



9k 



T k 

(63) 

where the measured electric field is given by E h sv . = 
E — V r Ep/e. In this situation we have, in addition to 



the electric current, a heat current (flux of heat per unit 
of area) given by 



U = — ^v k {e k -E F )g k 



(64) 



Both the electric and the heat currents can be written 
as21 

J = e 2 K ■ E obsv . + ^Ki ■ {-V r T) 

U = eK^E^. + ^Kz-i-VrT), (65) 

where Ki, i = 0,1,2 are second order tensors. In this 
problem the tensors are diagonal, i.e. Ki = lfcj, and by 
a well established procedure 2 ^ one obtains 



(\nE F /v ) 2 , 



(66) 



h = ~ (k B T) 2 ^ (InEp/io) 2 (l + (InEp/vo) x ) 



k 2 = ~(k B T) 2 ^-(\nEp/i ) 2 



(67) 
(68) 



In the above equations, we defined the two energy scales 
Uq = n riih v F and vq = Hvp/Rq. 

From the results (|66| , ((67|) and (|68|) it is easy to derive 
both the thermal conductivity k and the thermopower 
Q. These are given by 



1 r2 7r 2 



T 



3 h 



9 h 



(k B T) 2 ^-(\nEp/i Q ) 2 
u 2 

(k B T) 4 \ (1 + In E F /v Q ) 



and 



Q 



e 7 f T iE~ F ' 



(k B Tf (l + (In E F /to)- 1 ) 



(69) 



(70) 



Again, what should be emphasized in these results is 
the dependence of both k and Q on the particle den- 
sity, which is different from that of the usual two di- 
mensional electron gas and from the graphene sheet with 
only charged impurities in the substrate. Since it is ex- 
perimentally feasible to control the carrier density in the 
graphene plane 1 it is possible to check experimentally the 
dependence of the transport coefficients on the particle 
density. Finding the logarithmic corrections compared 
to the Coulomb scattering mechanism will be a strong 
indication for scattering due to midgap states. 

Normally, the second term of Eq. (f6"9"]l can safely be ne- 
glected and one obtains the the well-known Wiedemann- 
Franz law 



^ kfs T 



(71) 
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N = 2x10 4 









5x10 



I) 



5x10 



electronic density (cm ) 



electronic density (cm j 



FIG. 1: (color online). Conductivity as function of the elec- 
tronic density, for different values of temperature and for two 
impurity concentrations: (left) hi = 2 x 1CP 4 and (right) 
fn = 2 X 10" 3 . 
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at finite temperatures including scattering from charged 
impurities. Since there are two different mechanisms the 
total relaxation times is 



1 



D\k B T k _ Oi 
8h 2 pvgVF kp 



OL2kF , (72) 



where uo was defined in Eq. (|2U1) with Z = 1, e = 2.4 
and 7 = 4.2. 

The conductivity, including the contribution from two 
Dirac cones, reads 



e 2 E F k F 
X 4k B T 



x dx 



Q'l + Q^fcr 



• cosh 



2 ( EpX — fi 

2k B T 



where Ep = vphkp. The integral has a maximum around 
x = 1 and can be done numerically. The chemical poten- 
tial depends on the temperature and in the temperature 
range of T e [1,300] if is well described by the asymp- 
totic expression 



(irk B T) 2 
6E F 



(74) 



In Fig. [TJ the conductivity as function of the electronic 
density is shown for different values of temperature and 
for two impurity concentrations ra* = 2 x 10 -4 (left hand 
side) and ni = 2 x 10~ 3 (right hand side). In Fig. [2j 
the conductivity as function of temperature is shown, for 
two impurity concentrations fii — 2x 10~ 3 (upper panels) 
and fii = 2 x 10~ 4 (lower panels), and for two different 
electronic densities n — 7.2 x 10 12 cm~ 2 (left hand side) 
and n — 1.2 x 10 12 cm -2 (right hand side). 

For low impurity density, i.e., for realistic parameters, 
there is a striking temperature effect on the conductivity 
which is not seen in experiment. We, therefore, conclude 
once again that charged impurity can not resemble the 
main scattering mechanism. 



FIG. 2: Conductivity as function of temperature, for two 
impurity concentrations (up) hi — 2 x 10~ 3 and (down) 
fii = 2 x 1CP 4 , and for two different electronic densities (left) 
n = 7.2 x 10 12 cm" 2 and (right) n = 1.2 X 10 12 cm" 2 . 

But due to the logarithmic correction in the scatter- 
ing time, there is an additional term in ki, usually not 
present and thus a modified second term in Eq. (JSHJ). So 
even though our analysis is only valid for Ep/vq <C 1, 
we expect the Wiedemann-Franz law to be modified for 
large carrier densities. 

VI. NUMERICAL RESULTS 

A. Phonon contribution 

We now use the obtained values for the relaxation 
times of phonon scattering to compute the conductivity 



B. Influence of midgap states 

There are different estimates for the density of 
(charged) impurities. Preparing the Si-Si02 wafer by 
oxidizing the n-doped silicon-wafer produces no relevant 
impurity density which would affect the transport prop- 
erties of graphene. This is because charges due to dan- 
gling bonds are mainly localized between the Si-Si02 in- 
terface and thus exponentially suppressed due to the 300 
nm thick SiC>2 layer i 16 i 43 Placing the graphene sheet by 
micromechanical cleavage on top the wafer might pro- 
duce ionization of the OH-groups which neutralize the 
SiC>2-surface. We estimate a (relevant) charged impurity 
concentration of np < 10 11 cm _2 4i 

The estimates for the impurity density of vacancies is 
even lower having in mind the high energy cost of three 
missing bonds. Nevertheless, the main observation of this 
work is that midgap states give rise to a similar scatter- 
ing behavior as long-range Coulomb scatterers. Midgap 
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n=7.2x10 cm 




n j in [10 cm ] 



FIG. 3: (color online). The mobility fj, in units of cm 2 /Vs 
due to Coulomb scatterers with e = 3.9, 2.4; z — 5nm and 
e — 2.4; z — (solid lines), and vacancies as function of the 
impurity density for a carrier density n = 7.2 x 10 12 cm -2 
(dotted line). 



ity due to vacancies. Vacancies yield the lowest mobility 
for comparable impurity densities and thus represent the 
dominant scattering process. 

In Fig. 21 the conductivity is shown as function of the 
carrier density n. The left hand side shows the conduc- 
tivity due to Coulomb scattering with e = 2.4, z c = 
and nf — 2 x 10 10 cm~ 2 (dashed line), vacancies with 
m = 2 x 10 10 cm -2 (dotted line), and due to both con- 
tributions (full line). The right hand side shows the 
same plots with different Coulomb scatterers, i.e., with 
e = l,z c = 5nm and nf — 2 x 10 n cm~ 2 . 

The conductivity of the cleaner sample (left hand 
side) is not affected by the Coulomb scattering mech- 
anism and shows sub-linear behavior in the electronic 
density. For the dirty sample (right hand side) with 
z c = hnm, i.e., the charged impurities are well inside 
the substrate, the Coulomb scattering mechanism leads 
to super-linear behavior which in total yields linear be- 
havior of the combined conductivity with respect to the 
carrier density n. The above parameters yield mobilities 
of fi « 12000cm 2 /Vs (left) and fx « 8800cm 2 /Vs (right). 



-- Coul. (e=2.4, z=0, Nf=10 5 ) 


— Coul. (£=/, z=5nm, Nf=10~ 4 ) 


■■■■ MGS(N r 10 5 ) 
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FIG. 4: (color online). The left hand side shows the con- 
ductivity due to Coulomb scatterers with e = 2.4, z — 
and nf — 2 x 10 10 cm~ 2 (dashed line), vacancies with m = 
2 x 10 10 cm -2 (midgap states) (dotted line), and both contri- 
butions (full line). The right hand side shows the same plots 
with different Coulomb scatterers with e = 1, z — 5nm and 
nf = 2 x lO^cm -2 . 



states can also occur from cracks or boundaries. Another 
realization comes from impurities with large potential dif- 
ference with respect to the graphene sheet or corrugated 
graphene. We summarize all these effects in the impurity 
density n, < 10 11 cm -2 and set i?o ~ 1.4A. 

Fig. [3]shows the mobility /i = cr/ne at a carrier density 
of 7.2 x 10 12 cm -2 , which corresponds to a gate voltage 
of V g — 100V. The upper curves (solid lines) show the 
mobility due to Coulomb scattering, where the dielectric 
constant e and the average distance to the graphene sheet 
is varied. The lower curve (dotted line) shows the mobil- 



VII. CONCLUSIONS 

In summary, we presented a phenomenological the- 
ory for transport in graphene based on the semi- 
classical Boltzmann theory, first proposed by Nomura 
and MacDonald.® We pointed out that local point de- 
fects in form of vacancies, cracks etc. yield a similar k- 
dependence of the relaxation time as long-range Coulomb 
potentials. Moreover, they lead to a finite density of 
states at the Dirac point which can account for the ob- 
served minimal conductivity. The scattering mechanism 
due to midgap states has been widely ignored so far, but 
actually represents the dominant contribution to the to- 
tal conductivity. 

For "dirty" samples, this scattering mechanism yields 
an universal minimal conductivity, whereas for "cleaner" 
samples, the minimal conductivity increases. It also leads 
to a sub-linear behavior with respect to the carrier den- 
sity. In combination with Coulomb scattering, this be- 
havior may become linear. 

Regarding the numerical values, the major uncertainty 
lies in the impurities densities of charged or neutral de- 
fects. Cleaning the SiC>2 surface in a hydroxyl bat h 44 i 45 
would reduce charged impurities close to the graphene- 
sheet and thus estimate their effect to the conductivity. 
Another way of reducing a possible source of impurities is 
by interchanging the mechanical cleavage by "printing" 
the graphene sheet on top of the substrate 4£ 
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